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1 Introduction 

We summarize our work in collaboration with Maciej Zworski |lfi| . on the 
semiclassical density of resonances for a quantum open system, in the case 
when the associated classical dynamics is uniformly hyperbolic, and the set 
of trapped trajectories is fractal repeller. The system we consider is not a 
Hamiltonian flow, but rather a "symplectic map with a hole" on a compact 
phase space (the 2-torus). Such a map can be considered as a model for the 
Poincare section associated with a scattering Hamiltonian on M 2 , at some 
positive energy; the "hole" represents the points which never return to the 
Poincare section, that is, which are scattered to infinity. We then quantize 
this open map, obtaining a sequence of subunitary operators, the eigenvalues 
of which are interpreted as resonances. 

We are especially interested in the asymptotic density of "long-living reso- 
nances" , representing metastable states which decay in a time bounded away 
from zero (as opposed to "short resonances" , associated with states decaying 
instantaneously). Our results (both numerical and analytical) support the 
conjectured fractal Weyl law, according to which the number of long-living 
resonances scales as h~ d , where d is the (partial) fractal dimension of the 
trapped set. 

1.1 Generalities on resonances 

A Hamiltonian dynamical system (say, H(q,p) = p 2 + V(q) on R 2 ™) is said 
to be "closed" at the energy E when the energy surface Se is a compact 
subset of the phase space. The associated quantum operator H% = —H 2 A + 
V(q) then admits discrete spectrum near the energy E (for small enough K). 
Furthermore, if E is nondegenerate (meaning that the flow of H has no fixed 
point on Se), then the semiclassical density of eigenvalues is given by the 
celebrated Weyl's law |11|: 




(1) 
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This formula connects the density of quantum eigenvalues with the geometry 
of the classical energy surface E E . It shows that the number of resonances 
in an interval of type [E + Ch, E — Ch] is of order 0(h 1 ~"). Intuitively this 
Weyl law means that one quantum state is associated with each phase space 
cell of volume (2wh) n . 

When Ee is non-compact, or even of infinite volume, the spectral prop- 
erties of Hfi are different. Consider the case of a scattering situation, when 
the potential V(q) is of compact support: for any E > 0, Ee is unbounded, 
and Hfi admits absolutely continuous spectrum on [0,oo). However, one can 
meromorphically continue the resolvent (z — H%) across the real axis from 
the upper half-plane into the lower half-plane. In general, this continuation 
will have discrete poles {zj = Ej — i"fj} with "widths" jj > 0, which are the 
resonances of H^. 

Physically, each resonance is associated with a metastable state: a (not 
square- integrable) solution of the Schrodinger equation at the energy Zj, 
which decays like e~ tlj ^ h when t — > +oo. In spectroscopy experiments, one 
measures the energy dependence of some scattering cross-section cr(E). Each 
resonance Zj imposes a Lorentzian component ^_ g r ^ 2+ ^ 2 on o~(E); a res- 
onance Zj will be detectable on the signal o~(E) only if its Lorentzian is 
well-separated from the ones associated with nearby resonances of compara- 
ble widths, therefore iff \Ej — Ej\ ^S> Jj. This condition of "well-separability" 
is NOT the one we will be interested in here. We will rather consider the 
order of magnitude of each resonance lifetime fr/jj, independently of the 
nearby ones, in the semiclassical regime: a resonant state will be "visible" , 
or "long- living" , if jj = 0{h). Our objective will be to count the number of 
resonances Zj in boxes of the type {\Ej — E\ < Ch , 7j < Ch}, or equivalently 
{\zj-E\<CK}. 

1.2 Trapped sets 

Since resonant states are "invariant up to rescaling" , it is natural to relate 
them, in the semiclassical spirit, to invariant structures of the classical dy- 
namics. For a scattering system, the set of points (of energy E) which don't 
escape to infinity (either in past future) is called the trapped set at energy E, 
and denoted by K{E). The textbook example of a radially-symmetric poten- 
tial shows that this set may be empty (if V(r) decreases monotonically from 
r = to r — > oo), or have the same dimension as E E (if V(r) has a maximum 
V(ro) > before decreasing as r — > oo). 

For n — 2 degrees of freedom, the geometry of the trapped set can be more 
complex. Let us consider the well-known example of 2-dimensional scattering 
by a set of non-overlapping disks [5J|S| (a similar model was studied in |2"rJll2]). 

When the scatterer is a single disk, the trapped set is obviously empty. 

The scattering by two disks admits a single trapped periodic orbit, 
bouncing back and forth between the disks. Since the evolution between 
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two bounces is "trivial" , it is convenient to represent the scattering sys- 
tem through the bounce map on the reduced phase space (position along 
the boundaries x velocity angle). This map is actually defined only on a 
fraction of this phase space, namely on those points which will bounce again 
at least once. For the 2-disk system, this map has a unique periodic point 
(of period 2), which is of hyperbolic nature due to the curvature of the disks. 
The trapped set K of the map ("reduced" trapped set) reduces to this pair 
of points; it lies at the intersection of the forward trapped set _T_ (points 
trapped as t — > +00) and the backward trapped set r+ (points trapped as 
t — ► —00). 

The addition of a third disk generates a complex bouncing dynamics, for 
which the trapped set is a fractal repeller [Jj]. We will explain in the next 
section how such a structure arises in the case of the open baker's map. As 
in the 2-disk case, the bounce map is uniformly hyperbolic; each forward 
trapped point x £ /"!_ admits a stable manifold W-(x) (and vice- versa for 
x S J + ). One can show that n is fractal along the unstable direction W+: 
-T_ n W + has a Hausdorff dimension < d < 1 which depends on the positions 
and sizes of the disks. Due to time-reversal symmetry, J + n W- has the same 
Hausdorff dimension. Finally, the reduced trapped set K = r+ fl -TL is a 
fractal of dimension 2d, which contains infinitely many periodic orbits. The 
unreduced trapped set K(E) C Se has one more dimension corresponding 
to the direction of the flow, so it is of dimension D = 2d + 1. 



1.3 Fractal Weyl law 

We now relate the geometry of the trapped set K{E), to the density of 
resonances of the quantized Hamiltonian Hh in boxes { \ z — E\ < Ch}. The 
following conjecture (which dates back at least to the work of Sjostrand |2"T)1 
relates this density with the "thickness" of the trapped set. 

Conjecture 1. Assume that the trapped set K(E) at energy E has dimension 
2dE + 1. Then, the density of resonances near E grows as follows in the 
semiclassical limit: 

w ^ n #{Res(g ft )n{z : \z-E\<rh}} h ^ 

Vr > 0, — i j-j^ '- — ► c E (r) , (2) 

for a certain "shape function" < ce(t) < 00. 

We were voluntarily rather vague on the concept of "dimension" (a fractal 
set can be characterized by many different dimensions). In the case of a 
closed system, K{E) has dimension 2n— 1, so we recover the Weyl law QJ. If 
K(E) consists in one unstable periodic orbit, the resonances form a (slightly 
deformed) rectangular lattice of sides cx H, so each h-box contains at most 
finitely many resonances |22j . 

For intermediate situations (0 < d E < n — 1), one has only been able 
to prove one half of the above estimate, namely the upper bound for this 
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resonance counting |21M27llTUll24) . The dimension appearing in these upper 
bounds is the Minkowski dimension defined by measuring e-neighborhoods 
of K(E). In the case we will study, this dimension is equal to the Hausdorff 
one. Some lower bounds for the resonance density have been obtained as well 
|2"5] . but are far below the conjectured estimate. 

Several numerical studies have attempted to confirm the above estimate 
for a variety of scattering Hamiltonians jl()U12lH^lll4j . but with rather incon- 
clusive results. Indeed, it is numerically demanding to compute resonances. 
One method is to "complex rotate" the original Hamiltonian into a non- 
Hermitian operator, the eigenvalues of which are the resonances. Another 
method uses the (approximate) relationship between, on one side, the reso- 
nance spectrum of H^, one the other side, the set of zeros of some semiclassical 
zeta function, which is computed from the knowledge of classical periodic or- 
bits |SJ • I n the case 01 the geodesic flow on a convex co-compact quotient 
of the Poincare disk (which has a fractal trapped set), the resonances of the 
Laplace operator are exactly given by the zeros of Selberg's zeta function. 
Even in that case, it has been difficult to check the asymptotic Weyl law , 
due to the necessity to reach sufficiently high values of the energy |l()j . 

1.4 Open maps 

Confronted with these difficulties to deal with open Hamiltonian systems, we 
decided to study semiclassical resonance distributions for toy models which 
have already proven efficient to modelize closed systems. In the above ex- 
ample of obstacle scattering, the bounce map emerged as a way to simplify 
the description of the classical dynamics. It acts on a reduced phase space, 
and gets rid of the "trivial" evolution between bounces. The exact quantum 
problem also reduces to analyzing an operator acting on wavefunctions on 
the disk boundaries, but this operator is infinite-dimensional, and extracting 
its resonances is not a simple task [3] |2] . 

Canonical maps on the 2-torus were often used to mimic closed Hamilto- 
nian systems; they can be quantized into unitary matrices, the eigenphases 
of which are to be compared with the eigenvalues e~ lE i/ h of the propagator 
e -iH H /h j- gee e _g jyj and references therein for a mathematical introduction 
on quantum maps). 

We therefore decided to construct a "toy bounce map" on T 2 , with dy- 
namics similar to the original bounce map, and which can be easily quantized 
into an N x N subunitary matrix (where N = (2irh)^ 1 ). This matrix is then 
easily diagonalized, and its subunitary eigenvalues should be compared 
with the set {e~ 1Zj / h }, where the Zj are the resonances of Hn near some pos- 
itive energy E. We cannot prove any direct correspondence between, on one 
side the eigenvalues of our quantized map, on the other side resonances of a 
bona fide scattering Hamiltonian. However, we expect a semiclassical prop- 
erty like the fractal Weyl law to be robust, in the sense that it should be 
shared by all types of "quantum models" . To support this claim, we notice 
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that the usual, "closed" Weyl law is already (trivially) satisfied by quantized 
maps: the number of eigenphases 8j on the unit circle (corresponding to an 
energy range AE = 2wH) is exactly N = (2ttH)~ 1 , which agrees with the Weyl 
law for n = 2 degrees of freedom. Testing Conjecture Q in the framework 
of quantum maps should therefore give a reliable hint on its validity for more 
realistic Hamiltonian systems. 

Schomerus and Tworzydlo recently studied the quantum spectrum of an 
open chaotic map on the torus, namely the open kicked rotator |2()| : they 
obtain a good agreement with a fractal Weyl law for the resonances (despite 
the fact that the geometry of the trapped set is not completely understood 
for that map) . The authors also provide a heuristic argument to explain this 
Weyl law. We believe that this argument, upon some technical improvement, 
could yield a rigorous proof of the upper bound for the fractal Weyl law in 
case of maps. 

We preferred to investigate that problem using one of the best understood 
chaotic maps on T 2 , namely the baker's map. 

2 The open baker's map and its quantization 
2.1 Classical closed baker 

The (closed) baker's map is one of the simplest examples of uniformly hyper- 
bolic, strongly chaotic systems (it is a perfect model of Smale's horseshoe). 
The "3-baker's map" B on T 2 = [0, 1) x [0, 1) is defined as follows: 



This map preserves the symplectic form dq A dp on T 2 , and is invertible. 
Compared with a generic Anosov map, it has the particularity to be linear 
by parts, and its linearized dynamics (well-defined away from its lines of 
discontinuities) is independent of the point i e T 2 . As a consequence, the 
stretching exponent is constant on T 2 , as well as the unstable/stable direc- 
tions (horizontal/ vertical). 

This map admits a very simple Markov partition, made of the three 
vertical rectangles Rj = {q £ + l)/3), p £ [0,1)}, j = 0,1,2 (see 

Fig-D- Any bi-infinite sequence of symbols . . . e_2£-i ■ £o e i £ 2 ■ • ■ (where each 
e, £ {0, 1, 2}) will be associated with the unique point x s.t. B t (x) £ R et for 
all t £ Z. This is the point of coordinates (q,p), where q and p admit the 
ternary decompositions 




(3<z, §) if < g < 1/3, 

(3 ? -l,2+i)if l/3<g<2/3, 
(3g-2,4 2 )if 2/3<g<l. 



(3) 
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The baker's map B simply acts as a shift on this symbolic sequence: 

B(x = ... e_ 2 e-i • £o£ie2 ...) = ••■ ^-2^-1^0 ' ^2 ■■■ ■ (4) 

2.2 Opening the classical map 

We explained above that the bounce maps associated with the 2- or 3-disk 
systems were defined only on parts of the reduced phase space, namely on 
those points which bounce at least one more time. The remaining points, 
which escape to infinity right after the bounce, have no image through the 
map. 

Hence, to open our baker's map B, we just decide to restrict it on a subset 
S C T 2 , or equivalently we send points in T 2 \ S to infinity. We obtain an 
Anosov map "with a hole" , a class of dynamical systems recently studied 
in the literature 0]. The study is simpler when the hole corresponds to a 
Markov rectangle [3] , so this is the choice we will make (we expect the fractal 
Weyl law to hold for an arbitrary hole as well) . Let us choose for the hole the 
second Markov rectangle so that S = Rq U R%. Our open map C = B\s 
reads (see Fig. 0): 



if q £ Rq, 
\ ( 3g Z 2 , Z±2) if q e R 2 . 



(5) 



This map is canonical on S, and its inverse C 1 is defined on the set C(S). 




Fig. 1. Open baker's map C. The points in the middle rectangle are sent to infinity. 



Our choice for S coincides with the points x — {q,p) satisfying cq(x) <E {0, 2} 
(equivalently, points s.t. €q{x) — 1 are sent to infinity through C). This allows 
us to characterize the trapped sets very easily: 

the forward trapped set JL (see fig. EJ) is made of the points x which 
will never fall in the strip i?i for times t > 0: these are the points s.t. 
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ej G {0, 2} for all i > 0, with no constraint on the £j for z < 0. This set 
is of the form _T_ = Can x [0, 1), where Can is the standard 1/3-Cantor 
set on the unit interval. As a result, the intersection _T_ n W + = Can has 
the Hausdorff (or Minkowski) dimension d = ^| . 

the backward trapped set r + is made of the points satisfying 6 {0, 2} 
for all i < 0, and is given by [0, 1) x Can. 
the full trapped set K = Can x Can. 



Fig. 2. Iterative construction of the forward trapped set F- for the open baker's 
map C: we remove from T 2 the points leaving to infinity at times t = 1, 2, 3, 4 
(from light grey to dark grey) etc. At the end, there remains the fractal set _T_ = 
Can x [0, 1). 



2.3 Quantum baker's map 



We now describe in some detail the quantization of the above maps. We 
recall \7\ that a nontrivial quantum Hilbert space can be associated with 
the phase space T 2 only for discrete values of Planck's constant, namely 
h = (2ttN)^ 1 , N G No- In that case (the only one we will consider), this space 
Hn is of dimension N. It admits the "position" basis {Qj, j — 0, . . . , N — 1} 
made of the "Dirac combs" 



v). 



This basis is connected to the "momentum" basis {Pk, k 
through the discrete Fourier transform: 



(Pk\Qj) — {J~N)kj — 



^~2'nrNkj 



N 



j,ke{o,...,N-i}, 



,N- 1} 



(6) 
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where the Fourier matrix is unitary. Balazs and Voros proposed to 
quantize the closed baker's map B as follows, when N is a multiple of 3 (a 
condition we will always assume): in the position basis, it takes the block 
form 



This matrix is obviously unitary, and exactly satisfies the Van Vleck formula 
(the semiclassical expression for a quantum propagator, in terms of the clas- 
sical generating function). In the semiclassical limit N — > oo, it was shown |S] 
that these matrices classically propagate Gaussian coherent states supported 
far enough from the lines of discontinuities. As usual, discontinuities of the 
classical dynamics induce diffraction effects at the quantum level, which have 
been partially analyzed for the baker's map |25j (in particular, diffractive 
orbits have to be taken into account in the Gutzwiller formula for tr(_B^)). 
We believe that these diffractive effects should only induce lower-order cor- 
rections to the Weyl law 

We are now ready to quantize our open baker's map C of JSJ: since the 
classical map sends points in R\ to infinity and acts through B on S = 
Rq U i?2, the quantum propagator should kill states microsupported on 
and act as Bn on states microsupported on S. Therefore, in the position 
basis we get the subunitary matrix 



A very similar open quantum baker was constructed in |18j . as a quantiza- 
tion of Smale's horseshoe. In Figure El (left) we represent the moduli of the 
matrix elements (Cjv)nm- The largest elements are situated along the "tilted 
diagonals" n = 3m, n = 3(m — 2N/3), which correspond to the projection on 
the g-axis of the graph of C. Away from these "diagonals" , the amplitudes of 
the elements decrease relatively slowly (namely, like \/\n— 3m|). This slow 
decrease is due to the diffraction effects associated with the discontinuities of 
the map. 

2.4 Resonances of the open baker's map 

We numerically diagonalized the matrices Cat, for larger and larger Planck's 
constants N. First of all, we notice that the subspace Span{Qj , j = 
N/3, . . . ,2N/3 — 1}, made of position states in the "hole", is in the ker- 
nel of Cm- Therefore, it is sufficient to diagonalize the matrix obtained by 
removing the corresponding lines and columns. Upon a slight modification 
of the quantization procedure one obtains for CV a matrix covariant 
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3-baker (BV) 3-baker (W) 



0.5 m/N 1 0.5 m/N 1 




N= 27 N= 27 



Fig. 3. Graphical representation of the matrices Cat © and Cn 1101 . Each grey 
square represents the modulus of a matrix element (white=0, black=l/\/3) 

w.r.to parity, allowing for a separation of the even and odd eigenstates, and 
therefore reducing the dimension of each part by 2. This is the quantization 
we used for our numerics: we only plot the even-parity resonances (the dis- 
tribution of the odd-dimensional ones is very similar). In figure 0] we show 




-1 * ' 1 — ^ 1 1 ' ' 1 -1 1 1 ' ' I" — 1 1 1 1 

-1 -0.8 -0.6 -0.4 -0.2 0.2 0.4 0.6 0.8 1 -1 -0.8 -0.6 -0.4 -0.2 0.2 0.4 0.6 0.8 1 



Fig. 4. Even-parity spectrum of the matrices Cn for N/3 = 81, 2187 

the even-parity spectra of the matrix Cn for N — 3 5 and N — 3 8 . Although 
we could not detect exact null states for the reduced matrix, many among 
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the N/3 eigenvalues had very small moduli: for large values of N, the spec- 
trum of Cpf accumulates near the origin. This accumulation is an obvious 
consequence of the fractal Weyl law we want to test: 

Conjecture 2. For any radius 1 > r > and N £ No, 3\N, let us denote 

n{N, r) = #{A g Spec(CW) n { |A| > r } } . 
In the semiclassical limit, this counting function behaves as 

^ N -c(r), (9) 
N^ 

with a "shape function" < c(r) < oo. 

To test this conjecture, we proceed in two ways: 

• In a first step, we select some discrete values for r, and plot n(N, r) for 
an arbitrary sequence of N, in a log-log plot (see Fig. EJ . We observe that the 



Even eigenvalues of open 3-baker 

Counting eigenvalues in logarithmic scale 
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Fig. 5. Checking the JV-dependence of n(N, r) for various values of r, along geo- 
metric and arbitrary sequences for N. The thick curve has the slope log 2/ log 3 



slope of the data nicely converges towards the theoretical one (thick line) 
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Fig. 6. On the left, we plot the number n(N, r) of even eigenvalues of Cat of 

log 2 

modulus > r. On the right plot, we rescale those functions by the factors N lo s 3 



all the more so along geometric subsequences N — 3 fe N , and for relatively 
large values of the radius (r = 0.5). For the smaller value r = 0.03, the 
annulus {\z\ > r} still contains "too many resonances" and the asymptotic 
regime is not yet reached. 

log 2 

• In a second step, confident that n(N,r) scales like iV lo s 3 , we try to 
extract the shape function c(r). For an arbitrary sequence of values of N, we 
plot the function n(N, r) (Fig.EI left), and then rescale the vertical coordinate 

_ log 2 

by a factor N lo g 3 (right). The rescaled curves do roughly superpose on 
one another, supporting the conjecture. However, there remains relatively 
large fluctuations, even for large values of N. The curves corresponding to a 
geometric sequence N = 3 fc N Q , k = 0, 1, . . . tend to be nicely superposed to 
one another, but slightly differ from one sequence to another. Similar plots 
were given in 20- in the case of the kicked rotator; the shape function c(r) 
is conjectured there to correspond to some ensemble of random subunitary 
matrices. Our data are too unprecise to perform such a check. 

The fact that the spectra of the matrices Cn "behave nicely" along geo- 
metric sequences, while they fluctuate more strongly between successive val- 
ues of N, is not totally unexpected (similar phenomena had been noticed for 
the quantizations Bn of the closed baker P). In view of Fig. HO our conjecture 
(0 may be too strict if we apply it to a general sequence of N. At least, it 
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seems to be satisfied along geometric sequences { 3 fe N Q , k G N } , with shape 
functions cjq o {r) slightly depending on the sequence. 



3 A solvable toy model for the quantum baker 
3.1 Description of the toy model 

In an attempt to get some analytical grip on the resonances, we tried to 
simplify the quantum matrix Cat, keeping only its "backbone" along the 
tilted diagonals and removing the off-diagonal components. We obtained the 
"toy-of-the-toy model" given by the following matrices (the moduli of the 
components are shown on right plot of Fig- EJ : 



/l 1 








\ 


lOOOOOw 2 










lOOOOOw 










10 


1 







10 


J 1 







10 









10 





1 




10 





J 1 




\0 1 







) 



(10) 



From this example, it is pretty clear how one constructs Cat for N an arbitrary 
multiple of 3. A similar quantization of the closed 2-baker was introduced in 

Before describing the spectra of these matrices, we describe their propa- 
gation properties. Removing the "off-diagonal" elements, we have eliminated 
the effects of diffraction due to the discontinuities of C . However, this elim- 
ination is so abrupt that it modifies the semiclassical transport. Indeed, a 
coherent state situated at a point x away from the discontinuities will not 
be transformed by Cat into a single coherent state (as does Cat), but rather 
into a linear combination of 3 coherent states, shifted vertically by 1/3 from 
one another. Therefore, the matrices Cat do not quantize the open baker C 
of |J5J, but rather the following multivalued ("ray-splitting") map: 



C(a p)= /(39,!)U(3g,H±i)u(3 g ,H±!) 

\ (3q - 2, |) U (3g-2,£±i)U(3g -2, 



if q G Rq 



*±2) if qeR 2 . 



(11) 



This modification of the classical dynamics is rather annoying. Still, the dy- 
namics C shares some common features with that of C: the forward trapped 
set for C is the same as for C, that is the set n described in Fig. [3 On the 
other hand, the backward trapped set is now the full torus T 2 . 
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3.2 Interpretation of Cjv as a Walsh-quantized baker 

A possible way to avoid this modified classical dynamics is to interpret Cjv 
as a "Walsh-quantized map" (this interpretation makes sense when N = 3 fc , 
k € N). To introduce this Walsh formalism, let us first write the Hilbert space 
as a tensor product Hn = (C 3 )® fc , where we take the ternary decomposition 
of discrete positions jj = • eo£i • • • tk-i into account. If we call {eo, ei, 
the canonical basis of C 3 , each position state Qj E Hn can be represented 
as the tensor product state 

Qj = e eo ® e ex ® • • • ® e ek _ 1 . 

In the language of quantum computing, each tensor factor C 3 is the Hilbert 
space of a "qutrit" associated with a certain scale |19j . 

The Walsh Fourier transform is a modification of the discrete Fourier 
transform @, which first appeared in signal theory, and has been recently 
used as a toy model for harmonic analysis |15j . Its major advantage is the pos- 
sibility to construct states compactly supported in both position and "Walsh 
momentum" . In our finite-dimensional framework, this Walsh transform is 
the matrix 

(W N ) jr = 3 - fc / 2 exp(-^ ]T e*(Q£)), j,f = 0, . . . , N-l , 

t+i'=k-l 

and acts as follows on tensor product states: 

Wn {vq®vi®- ■ -Vk-i) = ^3Wfc-i®- • • FavigiFaVQ , ve € C 3 , £ — 0, . . . , fc-1 . 

Now, in the case N = 3 fc , our toy model Cjv can be expressed as 

fW N/3 \ 
c N = W^l o 

V w N/3 J 

One can show that "Walsh coherent states" are propagated through Cn 
according to the map C. Hence, as opposed to what happens in "standard" 
quantum mechanics, Cn Walsh-quantizes the open baker C. 

3.3 Resonances of C^ =3 k 

We now use the very peculiar properties of the matrices C 3 k to analytically 
compute their spectra. From the expressions in last section, one can see that 
the toy model Cjv acts very simply on tensor product states: 

Cn v <8> vi <g> • • • Vk-i = vi ® • • • v k -i ® J 7 3 ~ 1 tto2Vo , (12) 

where tt 2 projects C 3 orthogonally onto Spanjeo, e2}. Like its classical 
counterpart, Cn realizes a symbolic shift between the different scales. It also 
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Fig. 7. Spectra of the matrices Cn for N = 3 10 (circles), N = 3 15 (crosses). The 
large circles have radii |A| = 1, |A+|, A + A_|2. A_ 

sends the first symbol eo to the "end of the queue" , after a projection and 
a Fourier transform. The projection ttq2 kills the states Qj localized in the 
rectangle R\. The vector T% e eo in the last qutrit induces a localization in 
the momentum direction, near the momentum p = • eo. 

By iterating this expression k times, we see that the operator (C,/v) fe acts 
independently on each tensor factor C 3 , through the matrix J-^ 1 ttq2- The 
latter has three eigenvalues: 

it kills the state e\, implying that (Cjv) fe kills any state Qj for which at 
least one of the symbols ti{Qj) is equal to 1. These 3 fc — 2 k position states 
are localized "outside" of the trapped set _T_, which explains why they 
are killed by the dynamics. 

its two remaining eigenvalues A± have moduli |A + | w 0.8443, |A_| « 
0.6838. They build up the ^-dimensional) nontrivial spectrum of C/v, 
which has the form of a "lattice" (see Fig.[7|): 

Proposition 1. For N — 3 k , the nonzero spectrum of Cn is the set 

{A+}U{A_}|j{c 2i7r i A^ p/fe A^ /fc : 1 < p < k - 1 , 0<j<fc-l} . 

Most of these eigenvalues are highly degenerate (they span a subspace of di- 
mension 2 k ). When k — > oo, the highest degeneracies occur when p/k i=s 2. 
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which results in the following asymptotic distribution: 

V/eC(R 2 ), lim -j Yl mult(A) /(A) = / /" ( | A_ A+ 1 1 / 2 , 0) — . 

AeSpcc(C 3fc )\0 

The last formula shows that the spectrum of Cn along the geometric sequence 
{N — 3 fc , k S N} satisfies the fractal Weyl law ©, with a shape function 
in form of an abrupt step: c(r) = (5>(|A + A_ I 1 / 2 — r). Although the above 
spectrum seems very nongeneric (lattice structure, singular shape function), 
it is the first example (to our knowledge) of a quantum open system proven 
to satisfy the fractal Weyl law. 
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